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Abstract 

We develop a microscopic model that explains non-reciprocal optical effects 
in centrosymmetric Cr203. It is shown that light can couple directly to the 
antiferromagnetic order parameter. This coupling is mediated by the spin- 
orbit interaction and involves an interplay between the breaking of inversion 
symmetry due to the antiferromagnetic order parameter and the trigonal field 

contribution to the ligand field at the Cr^^ ion. 
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The study of the interaction of hght with magnetic substances has a long history. A 
classic example is the Faraday effect in ferromagnets where light couples directly to the 
ferromagnetic order parameter. Since the pioneering work of Argyres [jl|, it is known that 
electromagnetic radiation couples to the internal molecular field in a ferromagnet (which 
in turn is proportional to the order parameter) through the spin-orbit interaction. Such a 
coupling would of course be absent in antiferromagnets, where the internal molecular field is 
zero. In the absence of such a direct coupling between light and the antiferromagnetic order 
parameter, antiferromagnetic ordering could so far be probed only indirectly, for instance, 
by Raman scattering of magnetic excitations 0. 

The discovery of non-reciprocal optical effects (i.e., not invariant under time re- 
versal) 1^ below the Neel temperature, T/v, in optical experiments on Cr203 has therefore 
been considered a breakthrough in the study of antiferromagnetic ordering by light for it 
is only this class of experiments which can distinguish between two magnetic states that 
are related to each other by the time-reversal operation. Fiebig et al. IQ have found that 
antiferromagnetic domains could be observed directly by non-reciprocal second harmonic 
generation (SHG), leading to the first photographs ever of antiferromagnetic domains 
These experiments show that light can indeed couple directly to the antiferromagnetic order 
parameter. Though such a coupling was anticipated earlier from symmetry considerations 
P], no microscopic mechanism has been presented so far. 

In this Letter, we present a microscopic mechanism that explains all non-reciprocal op- 
tical effects in Cr203. While the spin-orbit interaction is, of course, essential in coupling 
the charge with the spin degrees of freedom it does not suffice for the generation of non- 
reciprocal effects. We find that non-reciprocal effects arise from an interplay between the 
breaking of crystal inversion symmetry by the antiferromagnetic order parameter and the 
trigonal distortion of the ligand field at the Cr^^ ion. This effect, in addition to the spin-orbit 
interaction, leads to a coupling of the antiferromagnetic order parameter with light. 

We present furthermore a simple cluster model, containing the full crystal symmetry 
of Cr203, which allows, for the first time, the orders of magnitude of all matrix elements 
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contributing to the non-reciprocal phenomena in Cr203 to be predicted. We apply the mi- 
croscopic model to the observed phenomenon of SHG and explain how antiferromagnetic 
domains can be distinguished experimentally. We also apply our model to another non- 
reciprocal effect seen experimentally in Cr203 viz., gyrotropic birefringence 0] and solve the 
long-standing question regarding its order of magnitude. 

As an introduction to SHG in Cr203, we discuss the macroscopic theory [Q in brief. 
Above Tn (— 307K), Cr203 crystallizes in the centrosymmetric point group D^d- The four 
Qj.3+ -Qj-^g ^Yie unit cell occupy equivalent c positions along the C3 (optic) axis. Since 
this structure has a center of inversion, parity considerations allow only magnetic dipole 
transitions related to the existence of an axial tensor of odd rank. Below T/v, time reversal 
symmetry (R) is broken and in this case, since Neumann's principle cannot be applied to 
non-static phenomena, only symmetry operations of the crystal that do not include R may 
be used to classify the allowed tensors for the susceptibilities. For Cr203, the remaining 
invariant subgroup is D-^. New tensors are allowed in this point group, for instance, a polar 
tensor of odd rank, that allow electric dipole transitions in SHG. 

From Maxwell's equations, one can derive the expression for SHG by considering the 
contributions to (non-linear) magnetization 

y^{2u,) ^ ^(2<^) . g(^) ^(uj) polarization 



p(2<^) 



X 



{2uj) 



■ Y^ii^) E*^'^). The source term S(r, t) in the wave equation [V x (Vx) -|- 



[1/c ) 9 /(9t ]E(r, t) = — S(r,t) can be written in a dipole expansion as |jTO 



, /92p(r,t) „ dM(r,t) 
S(r,t)=/io ( ^^ + Vx ^ + ... 



dt 



(1) 



Then, by assuming that E, P and M can be decomposed in a set of plane waves and 
considering a circular basis with E = -E+e^ + E^e^ + Eze^, where e± = =i=^(e^.±iey) and 
the direction of laser light to be along the optic axis, one obtains M, 



I 



-Im + ^Xe)E'i 
(7m + iXe)El 





(2) 
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where 7^ and Xe are non-zero components of the magnetic (7) and electric (x) susceptibihties 
that are allowed by D3 symmetry. Incoming right circularly polarized light (£'+) leads to 
left circularly polarized light and vice versa in SHG. Above T/v, as Xe = 0, the SHG 

intensities I± are identical while below T/v, the intensities I± oc | ±'-fm + iXe\'^E^ are different 
for right and left circularly polarized light, as observed in experiment 

The macroscopic theory dicussed so far is based on symmetry considerations only and 
does not provide estimates of the magnitude of either 7^ or Xe, which can be obtained only 
from a microscopic approach. A two ion mechanism for a non-zero electric dipole matrix 



element, Xe, in non-centrosymmetric antiferromagnets was proposed by Tanabe et al. []1T 
This mechanism is unlikely to be responsible for the coherent interference effect suggested 
by macroscopic theory (Q), since it involves magnetic excitations. Here, we propose instead, 
a single ion mechanism. In our theory the coupling of light to the antiferromagnetic order 
parameter is through the interference between the coherent contributions from the distinct 
Qj^3+ jQj^g crystal unit cell. No magnetic excitations are involved in our coupling 

mechanism. 

The crystal field due to the oxygen ions in Cr203 splits the five-fold degenerate 3(i- 
orbitals of the Cr^"*" ions into two levels, the lower one (^2 level) being triply degenerate 
{dxy, dyz, dzx) and the upper (e level) doubly degenerate (d^^-yi, d^^'^^r^). The three ^2 
orbitals are occupied in the ground state and the two e orbitals are empty. In SHG, the 
Qj^3+ jQj-^ absorbs two photons and is excited to a (^2)^6 configuration via two consecutive 
electric dipole processes (ED), corresponding to an r ■ E^"^) term in the Hamiltonian. A 
contribution to 7^ is then obtained for an emission via a magnetic dipole process (MD), 
corresponding to an L ■ B'^^'^^ term in the Hamiltonian. This contribution to 7m is allowed 
at all temperatures. The key point of a microscopic model is then to find the mechanism 
which allows a contribution to Xe via an electric dipole matrix element, r ■ E'-^'^^ in emission. 

We present our theory in two steps. As a first step, we explain the coupling mechanism 
between light and the order parameter and later, we proceed to discuss the role of the D^d 
crystal symmetry. In order to understand the origin of the ED transition below T/v, let us 



consider, for the moment, just two (i-orbitals per Cr^+ ion: the d^y orbital (ground state) 
and the dx2_y2 orbital (excited state). The MD contribution to SHG, i.e., to M*^^"^) is 

{d^y,ms\L\d^2_y2,ms){d^2_y2,ms\{v ■ E*^'^))^|(i^y, m^) , 

where the label denotes the spin quantum number of the relevant orbital. (Energy 
denominators have been omitted for clarity.) Here (r ■ E'-'^^)^ denotes r ■ E'^'^^ |q;) (a|r ■ E*^'^-' 
together with a sum over the intermediate states \q) and corresponding energy denominators. 

In order to develop the theory for the ED contribution, we treat the spin-orbit interaction 
and the trigonal distortion of the ligand field at the Cr^'^ ion as perturbations to the 3(i Cr 
states. It is easy to see that the diagonal part of the spin-orbit interaction mixes the t2 and 



e states, viz.. 



{d^y,ms\L-S\d^2_y2,ms) = 2i{Sz) , (3) 



where (S^) = (m^lS'^lm^). 



Next, we note that the 3d and the 4p Cr orbitals are mixed by the trigonal field [|12 
which breaks local parity. The local eigenstates of the Cr^"*" ion can therefore be written in 
the lowest order of perturbation theory as 

\dxy) = \dr,y,ms) + \{S^)\d^2_y2,ms) + T]'\p,ms) 
\dx2-y2) = \d^2_y2,ms) - \{S^)\d,^y,ms) +r]\p,ms) , 

where A is proportional to the spin-orbit coupling; rj, rj' are proportional to the trigonal field 
and \p,ms), \p',ms) are Cr Ap orbitals that will be specified later. The ED contribution to 
the SHG, i.e., to P^^"^), can now be expanded in powers of A and i] as 

{d^y\r\d^2_y2) {d^2_y2\{r ■ E^'^'>f\d^y) = 
r]X{d^y,ms\L ■ S\d^2_y2,ms){d^2_y2,ms\r\p,ms) (4) 
■(42_y2,m^|(r ■ E^'^^y\d^y,ms) + ... , 

where the contribution ~ Xi] is shown for illustration. All non zero-contributions are ~ 
Xri{Sz). To see this, consider all possible contributions order by order. 



The contribution ~ A"?]" vanishes because all d-orbitals have even parity, viz., 
{dxy,ms\r\dx2_y2,ms) = 0. The contribution ~ X'^r]^, viz., 

ri{d^y,ms\r\p,ms){d^2_y2,ms\{r ■ E^'^'>f\d^y,ms) (5) 

is finite for every Cr site and is proportional to rj. However, note that the trigonal field at Cr 
sites that are related by inversion symmetry have opposite signs, viz., if A1/B2 and B1/A2 
are pairs of Cr'^+ ions in the unit cell that are related by inversion symmetry, tja^ = —^B2 
and riA2 = ~VBi [Ol- summing the contributions to the matrix element ~ X^rj^ from 
each Cr ion in the unit cell, we see that this matrix element has to vanish identically. 

The term ~ X^rj^ vanishes, again because of the parity of the ci-orbitals but the contri- 
bution ~ X^rj^ shown in (|^) does not vanish generally. Using Eq. (^, we can sum up the 
total contribution ~ A^?7^ from the four (equivalent) Cr ions in the unit cell of Cr203 as 

-{S.)b2) = Ar^X^^-) A(T) , 

where X*^^"^) = {dx2-y2,ms\T\p,ms){d^2_y2,ms\{r ■ 'E^^^Y\dxy,ms) at site Al. Thus, we have 
shown that the ED contribution to the SHG, P^^"^), couples directly to the spin part of the 
antiferromagnetic order parameter A(T) = {Sz)ai — {Sz)bi + {Sz)a2 — {Sz)b2- Note that 
p(2t^) is non-reciprocal, as it changes sign under time-reversal and allows us to differentiate 
between the different antiferromagnetic domains in Cr203. 

We now complete the microscopic description by incorporating the full symmetry of 
Cr203, i.e., D^j^. This is done by considering a (Cr06)2 cluster model shown in Fig. |l]. In 
order to reproduce the full D^d symmetry, one has to choose the locations of the oxygen 
ions in this cluster in such a way that they do not coincide exactly with those in the actual 
crystal structure. This is because the cluster model contains only two Cr sites, while there 
are four in the unit cell of Cr203. 

We define the Cr sites to be located at Zo{0, 0, ±1) in the cluster model {zo is an arbitrary 
constant). Neglecting the trigonal distortion for the moment, let us choose the oxygens 
around the Cr sites to be located at 
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zo(0,0,e)+ ^6(7273,0,71/3) 

Zo(0, 0, e) + ^e(-7l76, /l/S) (6) 

Zo(0, 0, e) + ^9e(-yi76, -^1/2, v^l/S) , 

where = ±1 in order to obtain the sites of all the 12 oxygen ions, (see Fig. |l|). For any of 
the four combinations (e = ±1, -i? = ±1), the three oxygen sites given in (P) form a plane. It 
is obvious that this cluster model has a center of inversion at (0, 0, 0). In addition, it is easy 
to verify that this model has all the symmetry elements belonging to the group D^d 0, in 
particular, the C3 and the 2y symmetries that correspond to three- and two-fold rotations 
about the z- and y-axes respectively. 

The crystal field at the Cr site in the (Cr06)2 cluster is, in first approximation, cubic. 
The trigonal distortion of the ligand field can be treated as a perturbation. On doing this 
and rotating the crystal field axes with respect to the crystallographic axes (to facilitate 
calculations), we find that a convenient representation for the Cr one-particle orbitals is 
given by 

4^^ = -M3z2-r2) + Vl \Pz) 

4'^ = 3-V2[V2|42-j,2) + |4.)] + r/2b.) 

4'^ = 3-V2[^|4^) + |rf^^)] +r/2b,) (7) 

e(i) = 3-V2[|4.^)_V2|d,,)] + r/3b,) 

e(2) = 3-V2[|4,_^,)_y2|4.)] + r^sb,) , 
where the spin quantum numbers have been suppressed. For the first (second) Cr ion 
(e = ±1), one replaces z in the RHS of (|^) by {z =1= Zo) respectively. The above expression 
also includes the effect of the hemihedral part of the trigonal distortion which is the most 
dominant interaction. This interaction is of the form rjz and it leads to a mixing of Cr 3(i 
states and Cr 4p states with coefficients r/i, 772 and 773 being proportional to rj. Using (^, 
together with the spin-orbit interaction, one obtains after a lenghtly, but straightforward 
calculation, the dynamical current operator J(r, t) = (^! {t)\3\'^ {t)) . Here we report the 
results of this calculation, we will discuss the details elsewhere ||13|. The source term entering 
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is then obtained to be 



Since the full D^d symmetry of Cr203 is taken into account, the calculations based on the 
(Cr06)2 model correctly predict all selection rules also found by the macroscopic approach. 
In addition, we also obtain estimates for all non-zero matrix elements of the non-linear 
susceptibilities. 

The interference between the MD and the ED processes in SHG can be observed experi- 
mentally when the matrix elements, 7™, and Xe occuring in (|I]) are roughly of the same order 
of magnitude. We estimate their relative order of magnitude by 

— ~ 4 ^ TT-^ -eA^ A(T) . (8) 

Here Aq ~ 5000A is the wavelength of the emitted light, Oq ~ 0.69A is the radius of 
Cr^+, A ~ 100 cm~^ is the spin-orbit interaction [|14|,|15|, E^ — Et2 ~ SOOOcm"^ is the 



difference in energy between the ^2 and the e orbitals, rj ^ 350 cm ^ is the trigonal field 
\T^, Ep — Ed ^ 8 X 10^ cm~^ is the difference in energy between the d and the p orbitals 



that are mixed by the trigonal distortion and A(T) ^ 1 is the antiferromagnetic order 
parameter. The additional factor of 4 occurs since there are 4 matrix elements of the same 
order of magnitude (see @)). The above expression gives the right order of magnitude and 
we therefore conclude that the ED matrix element in this mechanism can indeed interfere 
with the MD matrix element. Clearly, the ED matrix element vanishes above T/v when time 
reversal symmetry is restored. 

We now consider the phenomenon of gyrotropic birefringence (GB). This is another non- 
reciprocal effect, the possible existence of which was first pointed out by Brown a/ 0. GB 
is a one-photon process appearing as a shift in the principal optic axes along with a change in 
the velocity of propagation of light. The first quantum mechanical treatment of this problem 
was presented by Hornreich and Shtrikman [|T^, who estimated that GB in Cr203 would lead 
to a shift in the optical axes of roughly 10~^ rad, viz., a very small effect. Recently however, 
Krichevtsov et al. measured this non-reciprocal rotation and the related magnetoelectric 
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susceptibility of Cr203 in the optical region 0. They found that the observed values were 
at least 4 orders of magnitude larger than those predicted by Hornreich and Shtrikman. 
They also found that the temperature dependence of the non-reciprocal effects mimicked 
that of the order parameter. The observed intensities and temperature dependence suggest 
that these effects originate from the ED process we have proposed. 

Now it is known that the dominant contribution to GB is from the magnetoelectric 
susceptibility defined by M*^"^^ = a : E*^"^) [T^. Using we have calculated the ED 
contribution to a in the optical region. We find that 

~ """^ "° ■ 

in dimensionless units. Here, Uo is the density of Cr ions in Cr203 (~ 3.3 x 10^^ m^^). In 
the region of experimental interest, h{uj — ujn) ~ 0.5 eV. Thus, we estimate a^x ~ 0.2 x 10~^ 
which is of the same order of magnitude as that observed experimentally. This also means 
that the non-reciprocal rotation would be ~ 10^^ rad. Since the ED process we consider 
couples light to the order parameter, the observed temperature dependence follows naturally 
from our mechanism. 

To conclude, we have developed a microscopic model that explains all non-reciprocal 
optical effects observed below Tjv in Cr203. We have shown that these effects can be ex- 
plained by an electric dipole process that arises from an interplay between the spin-orbit 
coupling and the trigonal distortion of the ligand field. Such a process couples light directly 
to the antiferromagnetic order parameter. Although we have applied the theory to explain 
non-reciprocal optical effects in Cr203, it can be generalized to all materials where i) the 
magnetic ion is not at a center of inversion and ii) inversion symmetry is broken below Tjv. 
In particular, we predict that such effects should be observed in V2O3, MnTi03 as also in the 
cuprate Gd2Cu04 below the ordering temperature of the Gadolinium magnetic subsystem. 



T7v(Gd)=6.5K p 
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FIGURES 

FIG. 1. An illustration of the (Cr06)2 cluster model. The circles and the triangles indicate the 
positions of the Oxygen and Chromium ions respectively. The direction of the triangles indicate 
the ordering of the Cr^+ moments in the antiferromagnetic state at a given domain. The cross is 
the location of the center of inversion. 
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